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Abstract 

Analysis of the average binary error probabilities (ABEP) and average capacity (AC) of wireless communica- 
tions systems over generalized fading channels have been considered separately in the past. This paper introduces 
a novel moment generating function (MGF)-based unified expression for the ABEP and AC of single and multiple 
link communication with maximal ratio combining. In addition, this paper proposes the hyper-Fox's H fading model 
as a unified fading distribution of a majority of the well-known generalized fading models. As such, we offer a 
generic unified performance expression that can be easily calculated and that is applicable to a wide variety of 
fading scenarios. The mathematical formalism is illustrated with some selected numerical examples that validate 
the correctness of our newly derived results. 

Index Terms 

Unified performance expression, average bit error rate, average capacity, maximal ratio combining, hyper-Fox's 
H fading channels, generalized Gamma fading, composite fading channels, extended generaUzed-K fading, moment 
generating function. 

I. Introduction 

The average binary error probabilities (ABEP) and average capacity (AC) are important performance 
metrics of wireless communication systems operating over fading channels. As such, considerable efforts 



IEEE TRANSACTIONS ON COMMUNICATIONS, VOL. X, NO. XX, JAN. 2011 



2 



have been devoted so far to develop analytical tools/frameworks to evaluate these two performance metrics 
and the references therein]. However, and to the best authors' knowledge, these tools/frameworks 
were developed separately and the computation of these two performance metrics was viewed as two 
independent problems. For instance, based on Craig's representation of the complementary error function, 
a unified moment generating function (MGF)-based approach was developed to compute the ABEP of 
a wide variety of modulation techniques over generalized fading |[T1 and the references therein]. More 
recently, other MGF-based approaches [|2l|-[|4l were also proposed for the capacity calculation of wireless 
channels subject to generalized fading. 

In contrast, this paper presents a novel MGF-based unified expression for the exact evaluation of both 
the ABEP and AC of single an multiple links generalized faded channels. The paper introduces also 
the hyper-Fox's H distribution as a versatile fading model including a variety of well-known models as 
special cases. With these two unifying frameworks at hand, we propose new generic expressions for the 
ABEP and AC with and without diversity reception. We also present some selected numerical examples 
to validate our newly derived results. 

The remainder of this paper is organized as follows. In Section [III a unified performance measure 
analysis of diversity receivers over additive white Gaussian noise (AWGN) channels is introduced and 
some key results are presented. In Section Uni new results for single link and multiple link reception 
are presented and applied to the newly proposed unifying hyper-Fox's H fading model. Finally, the main 
results are summarized and some conclusions are drawn in the last section. 

II. Unified Conditional Performance Expression 

A compact form for the conditional bit error probability (BEP) Pbep {lend) ^ certain value of 
instantaneous SNR 7^^^ for different binary modulations was proposed by Wojnar in [[5] Eq. (13)] as 



phase shift keying (PSK)), b depends on the type of detection for coherent, 1 for non-coherent), and 
r (-, ■) denotes the complementary incomplete Gamma function (61 Eq. (6.5.3)]. In the following theorem, 
we introduce an alternative representation of ^ using the incomplete beta function. 

Theorem 1 (Unified BEP Expression Using the Incomplete Beta Function). An alternative representation 




(1) 



where a depends on the type of modulation (| for orthogonal frequency shift keying (FSK), 1 for antipodal 
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of the compact form of the conditional bit error probability Pbep {lend) represented in ([l) is given by 
Ps.. = \ - ^^1^ f^jB (-^,,„,; M - </) . (2) 

where the parameters a and b depend on the particular form of the modulation and detection as mentioned 
above (see Tabled, and i = ^/—l is the imaginary unit with the property = —1, and B{z;a,b) = 
u'^~^{l — uf^^du is the incomplete beta function ^ Eq. (6.6.1)]. 

Proof: See Appendix |Al ■ 
In addition to the BEP performance measure, there exists another important performance measure 
commonly used in the literature, which is known as the conditional capacity. Explicitly, the conditional 
capacity is a measure of how much error-free information can be transmitted and received through 
the channel. The conditional normalized channel capacity Pc {lend) nats/s/Hz for a certain value of 
instantaneous SNR 7^^^ at the output of the receiver is well-known to be given by 

Pc {lend) = log (1 + a 7,„J , nats/s/Hz, (3) 

where a E M"*" represents the transmission power, and log {■) is, the natural logarithm (i.e., the logarithm 
to the base e) flU Eq. (4.1.1)]. We introduce in the following theorem a new alternative incomplete beta 
function-based representation of ([3]). 

Theorem 2 (Capacity Expression Using the Incomplete Beta Function). An alternative representation of 
the conditional capacity Pc {lend) given in dS]) is given by 

Pc{lend) = -B{-aJend-A,0). (4) 

Proof: See Appendix El ■ 
To the best of authors' knowledge, both ^ and ^ are not available in the literature. More importantly 
and using these two incomplete beta representations of the BEP and channel capacity measures, one can 
readily give, as shown in the following corollary, a unified performance expression whose special cases 
include the BEP and channel capacity. 

Corollary 1 (Unified Performance Expression Using Incomplete Beta Function). A compact and unified 
form of the conditional performance measure Pup {lend) (which include both the conditional BEP and 
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the conditional channel capacity) is given by 



Pup (lend) = « + P l^^j^^^ B [-^lend. hA - d) . (5) 

which reduces to © for a = 1, /3 = —1 and d ^ oo, and which reduces to ^ for a = 0, f3 = 2, a = 1, 
6=1 and d = 1. 

Proof: Based on Theorem \T\ and Theorem [2l the proof is obvious. ■ 
In the limit as d — oo, ([5]) might produce approximate results due to numerical computation limits 
of standard mathematical software packages. In this context, using both [T, Eq. (7.2.2/13)] and (Tl 
Eq. (7.3.1/28)], the unified expression given by Corollary [T] can be represented without limit as shown in 
the following corollary. 

Corollary 2 (Unified Performance Expression Using the Hypergeometric Function). The unified perfor- 
mance measure Pup {lend) '^^^ '^^^^ represented as 



" J 1 , / ,^n{a%nd) T? aW. 



Pup {lend) = 1 - 2 <) 1 + ("1) -^i [^"^ ^ + -«^e„d 



(6) 



which reduces to © for n = \, and which reduces to ([3]) for a = 1,6= 1 and n = 2. Moreover, in 



the coefficient set Ai™'' is defined as 



11 times 



Ai") = SXV^, (V) 



and [■; ■; ■] denotes the generalized hypergeometric function /0 Eq. (7.2.3/1)]. 

Proof: Using both [0 Eq. (7.2.2/13)] and El Eq. (7.3.1/28)], the proof is obvious. ■ 
Note that, in addition to the hypergeometric function representation given by Corollary [2l the unified 
expression can also be given in terms of other special functions such as the MacRobert's E function and 
the Meijer's G function, as shown in the following corollaries. These expressions are useful since they 
will facilitate in Section Un] the unified analysis of BEP and capacity in generalized fading environments. 

Corollary 3 (Unified Performance Expression Using the MacRobert's E Function). The unified perfor- 
mance measure Pup {lend) '^^^ ^^^^ be written as 

O'lend 

where E [■; ■; ■] is the MacRobert's E function ^ Sec. (9.4)], [|71 Sec. (2.23)]. Moreover, ^ reduces to 



(8) 
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([T]) for n = 1, and it reduces to ^ for a = 1 , 6 = 1 and n = 2. 

Proof: Using the generalized hypergeometric function representation of the MacRobert's E function 
[|7l Sec. (2.23)], i.e., substituting the following equality 

n^iTK) 



E [ai, a2, . . . , ar] bi, 62, • • • , bq] az] 



P Q 



Oi, 02, ... , ttr', 61, &2; • • • 5 bq] 



az 



(9) 



into © results in ([8]), which proves Corollary [3l 



Corollary 4 (Unified Performance Expression Using the Meijer's G Function). The unified performance 
measure Pup {lend) ^^■^ ^^^^ additional Meijer's G function-based representation given by 



Pup (7, 



n 



end) 



1) pl,n 



(^lend 



^(n) 

6,0 



where G™f 



is the Meijer's G function [|71 Eq. (8.3.22)]. Furthermore, (flOl) reduces to ©/or n 
and it also reduces to ^ for a = l,b = 1 and n = 2. 



(10) 

= 1, 



Proof: Using the relation between MacRobert's E and Meijer's G function [|71 Eq. (8.4.51/8)], we 
can write 



E 



1 



07, 



end - 



G 



l,n 
n,2 



^lend 



0,-6 



(11) 



Then, substituting (fTT)) into ^ results in (flOl) . which proves Corollary ID ■ 
It is worth to mention that the Meijer's G function is an special function defined by the Mellin- 
Bames type integral which contains products and quotients of the Euler gamma functions, and it is 
as such considered as a generalization of hypergeometric functions and other special functions such as 
exponential, bessel, logarithm, sine / cosine integral functions. Therefore, the PDF of some well-known 
fading distributions can be represented in terms of Meijer's G function. In this context, referring to 
Q Eq. (2.24.1/1)], the representation in Corollary |4] is more useful than the representations given in 
Corollaries [B |2] and [3] from numerical simplicity and computation start-points. 

In the following section, we use our new representations in the analysis of the BEP and capacity 
measures in generalized fading environments. 



III. Unified Performance Expression over Fading Channels 

Over the past four decades, the ABEP and the AC measures have been considered as two different 
problems, and many different solutions ranging from bounds to approximations, integral expressions, and 
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closed-form formulas have been presented for these two different measures for a variety of modulation 
schemes, diversity combining techniques, and fading distributions. At this point, we would also like to 
highlight again that these two different performance measures can be compactly combined as mentioned in 
the previous sections and can also be considered as a single problem. As such, for a certain nonnegative 
distribution of the instantaneous SNR 7^^^ at the output of the receiver (i.e., 7g„^ is distributed over 
(0, 00) according to P^/^^^i-f) which is the probability density function (PDF) of 7g„^), the average unified 
performance (AUP) expression is completely represented as 



where E[-] denotes the expectation operator, and as proposed in Section UIl Pup (7) is simply the 
conditional unified performance expression. 

A. Single Link Reception 

We consider an optimum receiver employing binary modulation and operating in a slow non-selective 
generalized fading environment corrupted by AWGN noise. In such a case, 7^^^ = 7^ is the instantaneous 
SNR at the output of the receiver, and it is distributed over (0, 00) according to the PDF p^^ (7). Referring 
to (fT2l) . the AUP expression can thus written as 



The PDF of the distribution of the fading is a nonnegative function (i.e., ^7^(7) > for 7 > 0) and a 
nonnegative function can be expressed in terms of Meijer's G or Fox's H function using integral transforms 
theory. In other words, the PDF of a variety of statistical envelope distributions, such as Rayleigh, Rician, 
exponential, Nakagami-m, WeibuU, generalized Nakagami-m, lognormal, K-distribution, generalized-K, 
etc., can be expressed in terms of either Meijer's G or Fox's H function [|9l. As a consequence of that, 
(fT3l) can be readily reduced to a closed-form solution by exploiting [|71 Eq. (2.24.1/1)]. In this context, 
the Fox's H channel fading model is introduced in the following definition. 

Definition 1 (Hyper-Fox's H Fading Channel). Let 7^ be a hyper-Fox's H fading distribution representing 
the instantaneous SNR 7^ at the output of the single-link receiver, which follows a probability law in 




(12) 




(13) 
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general such that its hyper-Fox's H PDF is given by 



{fnej, (f>na)j={ 



7>0, 



(14) 



with the conditions 



max 

i=l,2,..,M„^ 



< &n, < min 

i=l,2,..,Af„ 



nt e {l,2,...,Ki}, 



(15) 



such that (3 = ©1 n ©2 n . . . n &Ki, 7^ 0, where both {{fntj,<Pna)j2i} cind {{'&na-,Onij) j=i} coefficient 
sets are such parameters that they support the conditions given by [|7l Section 8.3.1], and where rjn^ G M 



and Cni, G M"*" are such two parameters that (14) certainly supports P'y^ il) = 1. In (fT4l) . H^^"[-] 
denotes the Fox's H function [jT] Eq. (8.3.1/l)lllc 

It is worth to notice that the PDF of most non-negative distributions can be compactly expressed or 
accurately approximated in the form of (fT4l) . For example, as seen in the first column of Tables IIIIIIIIVI 
and lYl the PDFs commonly used in the literature for modeling of instantaneous SNR distribution 7^^^ 
are listed in the form of hyper-Fox's fading distribution. From this perspective, where hyper-Fox's H 
distribution provides an unified framework on modeling of fading distribution, the unified and generalized 
result regarding the AUP expression of single-link reception can be obtained by means of substituting 
(fT4l) into ^ and using [El Theorem 2.9] with flVl Eq. (8.3.2/21)] such that the AUP expression over 
hyper-Fox's H fading channel is given by 



Tl 

Paup = 1 - - n 



;-ir 



M„^+n,Nni,+l 
Puf+'2,Quf+n 



^5 1); i'^-^riij "I" ^rifj; ^nij)j=\y (1)1; 



^mi))' 



(16) 



with the condition fl {©i fl ©2 fl . . . fl ©a-^} 7^ 0, where < ^ < 6 and the convergence region set 
{©„J is defined in (fTSl) . 

Let us consider some special cases of the hyper-Fox's H fading model in order to check analytical 
simplicity, accuracy and correctness of the AUP expression given by (fT3] ). 

Example 1 (Unified Performance Measure in Generalized Nakagami-m (GNM) Fading Channels). In 
GNM fading channels, the distribution of the instantaneous SNR 7^ follows a generalized Gamma PDF 

For more information about the Fox's H function, tfie readers are referred to fTOl. ifTTl 

^Using JT] Eq. (8.3.22)], tlie Fox's H function can be represented in terms of the Meijer's G function I?] Eq. (8.2.1)] which is a built-in 
function in the most popular mathematical software packages such as MATHEMATICA®. 
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given by 



Py, (7) 



r (m. 



(3, 



(17) 



for > 7 < 00, where the parameters > 1/2, > and 7^ > are the fading figure, the shaping 
parameter and the local mean power of the £th GNM distribution, and /S^ = F (m^ + 1/^^) /F (m^). It may 
be useful to notice that the special or limiting cases of the GNM distribution are well-known in the literature 
as the Rayleigh (m^ = 1,,^^ = 1), exponential (m^ = 1,,^^ = 1/2), Half-Normal (m^ = 1/2,,^^ = 1), 
Nakagami-m (^^ = 1), Gamma (^^ = 1/2), WeibuU (m^ = 1), lognormal (m^ — )■ 00,^^ — )■ 0), and 
AWGN (m^ — 00,^^ = 1). Referring to Definition [T] in order to obtain the AUP expression Paup for 
generalized Gamma fading channels by means of using (fT6l) . one can readily represents the PDF (7) 
of the generalized Gamma distribution in terms of the hyper-Fox's H distribution by means of using [[TTl 
Eq. (2.9.4)], i.e., using the Fox's H representation given in the first row of Table IVl 



(7) 



/3. 



1^ i'^i 



tt1,0 



'^7 
,7£ 



i/e„i/e,)J 



(18) 



where — means that the coefficients are absent. Next, expressing (1181) in terms of ([14)) with K = 1, and 
then substituting (fTSi) into (fT6l) (i.e., by means of mapping the parameters of two PDF model (fT4l) and 
(fT8l) . and then utilizing (fT6l)). one can readily obtain the AUP expression for GNM fading channels as 



TTn+1,1 
■n-2,n+l 



«7£ 



(1-^1), (1,1) 
A[o,\),(m„l/e,) 



(19) 



Here, it might be useful to notice that, in addition to this AUP expression for GNM fading channels. 
Table in offers simplified expressions for the AUP of a variety of commonly used fading channels in order 
to facilitate for the readers the use of our unified BEP and AC results. In order to check the validity and 
completeness of (fT9l ), substituting n = 1 results in the ABEP of signal transmission over GNM fading 
channels as expected, that is 



Pabep ~ '2^^ 



-H 



2,1 
2,2 



L«7£ 



(i-&,i),(i,i: 

(0,l),(m„l/e, 



(20) 



V{h)V{m, 

After performing some algebraic manipulations either by means of using [[TTl Property 2.11] or by 
recognizing that the cumulative distribution function (CDF) of a Fox's H distribution can be written 
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in two different forms given in lfT2l Eq. (4.17)] and |[l2l Eq. (4.19)], (|20l ) can be simplified to 



P, 



ABEP 



-H 



1,2 



(1,1), (1-^1) 

(m„l/e,),(0,l) 



(21) 



2T{h)T{m,y 

in perfect agreement with |[T3l Eq. (10)]. Note that in the special case of Nakagami-m fading channels, 
substituting = 1, it is straight forward to show that (l20l) reduces to [5, Eq. (17)] by means of some 
algebraic manipulations using [|7J Eq. (8. 4.49/13)] and [|71 Eq.(7.3.1/28)] together. In addition, for Rayleigh 



fading channel (i.e., = = 1), (1201) simplifies to Pabep 



which agrees with the 



2 I yi+a-y^^ 

four results given in [[D for BFSK (a = 1/2 and h = 1/2), BPSK (a = 1 and 6 = 1/2), non-coherent 
BFSK (NC-BFSK) (a = 1/2 and 6=1) and differentially encoded BPSK (BDPSK) (a = 1 and 6 = 1). 
In addition to these expected ABEP consequences of (fT9l) . the AC of the GNM fading channels can be 
also easily obtained by setting 6=1 and n = 2 in (fT9l ), yielding 



AC 



1 tt3,1 



«7£ 



(0,1), (1,1) 
(0,l),(0,l),(m„l/e,) 



(22) 



where a denotes the transmitted power. Setting = 1 and then using [|7l Eq.(8.3.2/21)], (l22l) further 



reduces to the well-known result P^c = G 



3,1 
2,3 





0,1 




0,0,m^ 



/r(m^) which is the AC of the Nakagami-m 



fading channels lfT4l Eq. (3)]. Furthermore, using IfTTl Eq. (2.1.7)] and [[7J Eq. (8.4.11.3)], then recalling 
the relation between the first order En integral Ei{x) and exponential integral Ei{x) such as Ei{x) = 
—Ei{—x), for Rayleigh fading channel (m^ = = 1), (l22l) simplifies to [HI Eq. (15.26)] as expected. 

In addition to the review above regarding the unified performance measures in GNM fading channels, 
let us consider the shadowing identified as a main cause not only for reducing energy but also causing 
performance loss and instability at the receiver. Specifically, shadowing is the effect that the received 
signal power, i.e., the mean of instantaneous SNR fluctuates due to objects obstructing the propagation 
path between transmitter and receiver. Additionally, in the architecture of next- generation wireless commu- 
nication systems, a spectrum well above 60 GHz frequency will be used, causing greater susceptibility to 
shadowing. In this context, the modeling of shadowing plays an important role in the context of designing 
systems and evaluating the corresponding performance. The most well-known shadowing distribution in 
the literature is the Lognormal distribution [[H and the references therein]. 

Example 2 (Unified Performance Measure in Lognormal Fading Channels). The PDF of lognormal 
distribution can be well approximated in the form of (fT6l) as shown in the second row of Table |IV] 
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such that 

I \ '"^fe ttO,0 

where /if(dB) and crf(dB) are the mean and the standard deviation of 7^, and where Uk is defined as 
ojk = lo('^''e^k+iJ.e)/io gy^,]^ jj^^j ^Qj. ^ ^ {1, 2, ... , K}, {wk} and {uk} are the weight factors and the zeros 
(abscissas) of the K-order Hermite polynomial Table 25.10]. Mapping the parameters and coefficients 
between (|23|) and (fT4)) . one can readily give the unified performance measure over lognormal fading 
channels by means of utilizing (fT6l) and then using [|7l Eq. (8. 3.2/21)] as 

K 

As expected, setting n = 1 in (|24)) results in Pabep = 2^ Yl ^fcr(&, aa;fc)/r(6). Moreover, setting 6 = 1 

71=1 

and 7T, = 2, (|24l) simplifies into Pabep = S ^fc ^^S + '^^k) as expected. 

71=1 

In addition to the context of shadowing, composite fading channels (include both shadowing and fading) 
plays an important role in designing and modeling wireless communication systems. To the best of our 
knowledge, another general composite fading model is the extended generalized-K (EGK) fading. In the 
following example, unified performance measure Paup is analyzed over EGK fading channels ifTSl . 

Example 3 (Unified Performance Measure in Extended Generalized-K (EGK) Fading Channels). The 
distribution of the instantaneous SNR in the EGK fading channel follows the PDF given by [[TSl Eq. (3)], 
lfT6]| . that is, 

defined over < 7 < 00, where the parameters (0.5 < < 00) and (0 < < 00) 
represent the fading figure (diversity severity / order) and the fading shaping factor, respectively, while 
f^si (0-5 < < 00) and (0 < < 00) represent the shadowing severity and the shadow- 
ing shaping factor (inhomogeneity), respectively. In addition, the parameters (3^ and (3^^ are defined as 
= r (m^ + 1/^^) /r (m^) and = T (m^^ + 1/^^) /F (m^^), respectively. In addition, F (-, -, -, ■) is 
the extended incomplete Gamma function defined as F (a, x, b, (3) = r"~^ exp (— r — 6r~'^) dr, where 
a, /3, 6 G C and x G M"*" [[TtI Eq. (6.2)]. The Fox's H representation of (|25]) is given by at the second row 



7 









(23) 



auk 



l-b,l 



A 



(n) 
(0) 



(24) 
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of Table |Vl that is 



Pi, (7) 



r(m^)r(mj7^ 



H 



2,0 
0,2 



7, 



'-1 



si 



(m^ - 1/^^, 1/^^), {m^^ - 1/4^, 1/^,^) _ 



(26) 



Mapping the parameters and coefficients of (l26l) to those of the hyper-Fox's H fading model defined in 
Definition [T] with K = \, the unified performance measure Paup can be readily expressed as 



P 



n 



'-11 



AUP 



TTn+2,1 



Mil 

^Isl 



(1-&,1),(1,1) 



^((0,1))' ^^i^ i 



(27) 



r(6)rK)rK 

Note that (1271 ) reduces to the ABEP in EGK fading channels [fT6l Eq. (34)] by means of setting n = 1 
and using the same steps in the derivation of (|2TI) . Moreover, for n = 2 and 6=1, (|27] ) simplifies into 
the AC in EGK fading channels [[T6l Eq. (42)] as expected. 

In order to check the analytical accuracy and correctness, some numerical and simulation results 
regarding the ABEP and AC performances of single link reception over fading channels are depicted 
in Fig. m Fig. |2] and Fig. |3l and they show that these numerical and simulation results are in perfect 
agreement. 



B. Multiple Link Reception 

We consider an L-branch maximal ratio combiner (MRC) diversity system employing binary modulation 
and operating in a slow non-selective mutual independent and not-necessarily identically distributed 
generalized fading environment corrupted by AWGN noise. The instantaneous SNR 7^^^ at the output of 
the MRC receiver is considered as the sum of the instantaneous SNRs of the branches, that is, 

L 

lend = Yl ^^^^ 

where L denotes the number of branches, and where for i E {1, 2, 3, ... , L}, 7^ is the instantaneous SNR 
the ith branch is subjected to. The AUP expression Paup of the L-branch MRC combiner can be obtained 
by averaging the (instantaneous) unified performance measure Pup {l^nd) given by Corollary |4]over the 
PDF of 7e„rf = X]^=i 1i shown in (fT2l) . Due to several reasons (e.g., insufficient antenna spacing or 
coupling among RF layers), correlation may exist among diversity branches of the L-branch MRC. With 
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that, the AC using (fT2l) involves an L-fold integral given by 



71 

Paup = 1 - 2 ^ ^ 



-1)' 



n poo poo 



m Jo Jo 



pl,n 



L 

1 



in) 



X 



L-fold 



6,0 



rL)dr^dr2. . .drL}, (29) 



where ^7^,72,.. .,7^ ('"i, ■ ■ ■ , r^) is the joint multivariate PDF of the instantaneous SNRs {7^}^^^ In (|29l) . 
the L-fold integration is tedious and cannot be partioned into the product of one dimensional integrals 
even if the instantaneous SNRs are assumed mutually independent. Additionally, it is clear that 

the numerical evaluation of (|29l ) is complex requiring a long time to compute the desired result as the 
number of branches L increases. Fortunately, after performing some algebraic manipulations, the AUP 
expression can be readily obtained in terms of a single integral expression using an MGF-based approach 
[[U as presented in the following theorem. 

Theorem 3 (Average Unified Expression of the L-branch Diversity Combiners over Correlated Not-Nec- 
essarily Identically Distributed Fading Channels). The exact AUP of L-branch diversity combiner over 
mutually not-necessarily independent nor identically distributed fading channels is given by 



P. 



AUP 



n 
2 



G 



l,n 

n+1,2 



in) 



6,0 



ds 



(30) 



where the parameters a G b G (0,1], n G {1,2} are selected according to desired performance 

lere M^^^^^^,„^^^ {s) i 
SNRs 7]^, 72, . . . , 7^ of the branches. 



measure, and where A^7j,72,...,7^ (s) = E[exp(— s 'Yl,ilt)\ joint MGF of the correlated instantaneous 



Proof: See Appendix O ■ 
It is worth accentuating that, in order to find the AUP expression of the diversity combiner, the MGF- 
based technique proposed in Theorem [3] eliminates the necessity of finding the PDF of the instantaneous 



SNR 7g„^ = J2e=i li through the inverse Laplace transform (lUT) of the joint MGF M 



7i,72v,7i 



E[exp(— s J2e1e)]- Shortly, Theorem [3] suggests that one can readily obtain the AUP expression using the 
joint MGF (s). Additionally, the integral in (|30l ) can be accurately estimated by employing 

the Gauss-Chebyshev Quadrature (GCQ) formula [H Eq.(25.4.39)], i.e.. 



AUP 



n 
2 



(-1)" 

m 



N 



n=l 



a 


aS"\i" 




Sn 


6,0 


{ 



d 

^-M7i,72,-,7^ (^) 



(31) 
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which converges rapidly and steadily, requiring only few terms for an accurate result. In (|3T| ). the 
coefficients Wr, and s„ are defined as 



TT sm 



Wr. 



\ 2N 'V 



and 



5„ = tan(fcos(^7r) + f), (32) 



4Arcos2 (|cos(^7r) +f) 

respectively, where the truncation index could be chosen more than = 30 to obtain a high level of 
accuracy. Despite the fact that the novel technique represented by Theorem [3] is easy to use, and referring 
to both (l30l) and (|3TI) . let us consider its special cases in order to check its analytical simplicity and 
accuracy. 

Special Case 1 (Average Bit Error Probabilities of L-Branch Diversity Combiner). As mentioned before, 
the ABEP Pabep = |lE[r(6, a7g„^)/r(6)] of the L-branch diversity combiner can be readily obtained 
through setting n = 1 in (|30l) . such that Pup\n=i results in 



ABEP 



G 



1,1 

2,2 



1,1 

6,0 



ds. 



(33) 



2 2T{h) J, 

To the best of the authors' knowledge, the result given by (l34l) is a new result which readily simplifies 
into the performance of BFSK (a = 1/2 and h = 1/2), BPSK (a = 1 and 6 = 1/2), non-coherent BFSK 
(a = 1/2 and 6=1) and BDPSK (a = 1 and 6 = 1). 

Special Case 2 (Average Capacity of L-Branch Diversity Combiner). The AC Pac = IE[log(l + a7g„^)] 
of the L-branch diversity combiner can be readily obtained through setting 6 = 1 and n = 2 in (l30l) . such 
that Pup\n=2b=i results in 



AC 



ni.2 

*-'3,2 



1,1,1 
1,0 



ds, 



(34) 



where it may be useful to notice that f ^[\^ = GJ^ [f | ^ ] by means of employing Eq. (8.2.2/9)] . 
Then, using both [7, Eq. (8.2.2/14)] and G2;i[f | Y] = -Ei (-f) E Eq. (8.4.11/1)], where Ei (•) is the 
exponential integral function ^ Eq. (8.211/1)], (|34l) simplifies to 

d 



AC 



Ei{-s/a) 



is) 



ds. 



(35) 



which is a well-known result given by Di Renzo et al. in [|3l Eq. (7)]. 



Note that the spatial correlation between all fading amplitudes can be determined from the physical 
parameters of the model, which includes antenna spacing, antenna arrangement, angle spread, and angle 
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of arrival. In the case of there does not exist any correlation between all fading amplitudes 7^^, 73, . . . , 7^ 
for the branches of the diversity combiner, the AUP expression is given in the following corollary. 

Corollary 5 (Average Unified Expression of the L-branch Diversity Combiners over Mutually Independent 
Non-Identically Distributed Fading Channels). The exact AUP expression of L-branch MRC over mutually 
independent and non-identically distributed fading channels is given by 



n poo 



AUP 




p l,n 
^n+1,2 



6,0 



=1 



l[M^^{s)ds\, (36) 



-> fc=i 



where, for i E {1, 2, ... , L}, -M-y^ (s) = E[exp (—57^)] is the MGF of the instantaneous SNR 7^ that the 
ith branch is subjected to. 

Proof: When there is no correlation between all instantaneous SNRs 7^^, 72, . . . , 7^, of the branches, 
one can readily write A^-y^,72v,7l('^) ~ Il^li -^7^ i^) whose derivation with respect to s is 



L r 



7fe 



(37) 



Finally, substituting (l37l) into (l30l) results in (l36l) . which proves Corollary [51 ■ 
As mentioned before and nicely shown in Tables [III |llll |IV] and |Vl the PDF of several non-negative 
distributions can be compactly expressed or accurately approximated in the form of (fT4l) . Using either 
lfT2l Eq. (2.12)] or lUl Eq. (2.10)], the MGF of the hyper-Fox's H fading channel is given by 



rim j^Mn^,Nni,+l 

n,=l 



'^nijjj=i 



(38) 



with the convergence region 3ft {s} > 0. Then, with the aid of [|7l Eq. (8.3.2/15)], the derivative of (|38] ) 
can be readily derived as 



d 



ng=l 



(0, 1), {^mj + ^na, Ornj)jJi 
i'-fnij + 4'ntj, 4>nij)j=i 



(39) 



Finally, substituting both (1381) and (1391) into (1361) . the AUP expression of the L-branch diversity combiner 
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performing over hyper-Fox's H fading channel can be readily obtained in the form of Corollary [5] as 



AUP 




l,n 
n+1,2 



b,0 



X 



EE 

£=l ne=l 



-rig 



(0,l),(^n,i 



k=i nfc=l 



'nijj (Pntj)j=l 



ds 



(40) 



Let us consider some special cases of the hyper-Fox's H fading model, i.e., the special cases of (l40l) . 
Note that the MGFs of some commonly used fading distributions (such as one-sided Gaussian, expo- 
nential. Gamma, WeibuU, hyper-Gamma, Nakagami^ (Hoyt), Nakagamin (Rice), Maxwell, lognormal, 
K-distribution, generalized- K, generalized Gamma, extended generalized Gamma and Fox's H) and their 
derivatives are given in details in Tables HHIIIIIVI and IVl Using the MGF of generalized Gamma distribution 
and its derivative given in the first row in Table IVl and then substituting them into (|40|) . the AUP expression 
of the L-branch combiner over GNM fading channels is given by 



AUP 




(-1)" 

Tib) 



n 

.k=l 



-G 

s 

EHu 



l,n 
n+1,2 



7^5 



Af'M 



(n) 

1 5 

6,0 
(0,1) 



X 



k = l 

k^e 



(1,1) 



ds 



(41) 



Note that in the special case of single link, i.e., L = 1, (|4TI) not surprisingly simplifies into (fT9l) by means 
of some algebraic manipulations using |l7l Eq. (2. 25. 1/1)] and [7, Eq. (8. 3. 2/21)] together. Some algebraic 
manipulations substituting the shaping parameters = ^2 — ' ' ' — — ^ ^^^^ (l4T]) and then utilizing 
[|7l Eqs. (8.3.2/7), (8.3.2/21) and then (8.4.2/5)] result in the unified expression of the L-branch diversity 
combiner performing over mutually independent and not necessarily identically distributed Nakagami-m 
fading channels, that is. 



n 



AUP 



G 



l,n 

n+1,2 



a!"\i 

6,0 



L 

E 



7£ 



1 + 



n 1 + ^ 



-ds 



(42) 



For identically distributed Nakagami-m fading channels (i.e., m-^ 



rrir, 



1l 



7), apparently as a result of both employing V{mL + 1) (l + ^) 



g \ —mL—l 



m and 7i = 72 = 

gl,l I —rriL 1 
1)1 L m I J 
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Id Eq.(8.4.2/5)] and then using the integral equality of two Meijer's G functions IH Eq.(2.24.1/1)], ^ 
simplifies to 



n 



m 



a7 



1 - 6, 1 



(43) 



which is the special case of (fT9l ) with = 1, = m, 7^ = 7 and L = 1 as it is expected. In addition, 
for n = 1, following the same steps in the derivation of (|2TI) from (l20l) . the unified expression given by 
(l42l) readily reduces to the ABEP of the L-branch diversity combiner over identical Nakagami-m fading 
channels can be readily obtained as 



ABEP 



-G 



1,2 



m 


1,1-6" 


07 


0, rriL 



2T{b)T{mL) 

For n = 2 and 6=1, (|43] ) simplifies to the well-known result flU Eq. (33)], that is. 



(44) 



1 



AC 



3,1 



r(mL) ^2'=^ 





0,1 


a7 


0,0, mL_ 



(45) 



where a denotes the transmitted power. 

As an illustration of the mathematical formalism presented above, some numerical and simulation results 
regarding the ABEP and AC performance of multiple link reception over fading channels are depicted in 
Fig. m Fig. [5] and Fig. [6l and these figures show that these analytical and simulation results are in perfect 
agreement. 



IV. Conclusion 

In this paper, we presented a unified performance expression combining the ABEP and AC of wireless 
communication systems over generalized fading channels. More precisely, this paper introduces an MGF- 
based unified expression for the ABEP and AC of single and multiple link communication with an L-branch 
MRC combining. In addition, the hyper-Fox's H fading model is proposed as a unified fading distribution 
for a majority of the well-known generalized fading models in order to provide more general and more 
generic results which can be readily simplified to some published results for some well-known fading 
distributions. We explicitly offer a generic unified performance expression that can be easily calculated 
and that is applicable to a wide variety of fading scenarios. Finally, as an illustration of the mathematical 
formalism, some simulations have been carried out for different scenarios of fading environment, and 
numerical and simulation results were shown to be in perfect agreement. 
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Appendix A 

Unified BEP Expression Using the Incomplete Beta Function 

Note that, substituting the lower incomplete gamma 7 (a, z) = e~^u"-~^du ^ Eq. (6.5.2)] and using 
the relation between the lower and upper incomplete Gamma functions (i.e., 7 (a, z) = T (a) — T (a, z)) 
Il6l Eq. (6.5.2)], the conditional bit error probability Peer (lend) '■^^^ ^e written as 

PbEP [lend) = 2 2r^) ■ ^ ' 

Substituting an alternative representation of the lower incomplete gamma function 7 (a, z) = h~^{az)^ ^F^ [a; a + 1; 
lIT] Eq. (7.11.3/1)] in ( lA.ll) and then using the well-known limit representation ^F^[—z]a]a + l\ = 
lim 2-^1 [d, a; a + 1; —z/d] E Eq. (7.2.2/13)], we can write 



^B.P (7„.) = i - llfj. Ita M H- l; -^7„J , (A.2) 



where ^F-^^ [■; ■; ■] and 2-^1 [s ") "i ■] the Kummer confluent and Gauss hypergeometric functions, re- 
spectively. Finally, substituting the Gauss hypergeometric representation of the incomplete beta function, 
i.e., 2^1 [z;d,b;b+ 1] = bz^''B {z;b,l - d) [0 Eq. (7.3.1/28)] into (IXIl) results in ©, which proves 
Theorem [U 

Appendix B 

Capacity Expression Using the Incomplete Beta Function 
Using (H Eq. (7.3.3/7)], i.e., 

n-l 

inn J n 

fc=0 

one can readily show by setting m = = 1 in (IB. II) that the conditional capacity Pc [lend) ^'^^ 
represented as 

Pc [lend) = lend 2^1 [1, 1; 2; "7^ J • (B.2) 



Substituting the Gauss hypergeometric representation of the incomplete beta function, i.e., 2-^1 [z; a,b;b + 1] = 
bz~^B [z] 6, 1 - a) l|7l Eq. (7.3.1/28)] into (lR2l) results in @i, which proves Theorem H 
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Appendix C 

Average Unified Expression of the L-branch Diversity Combiners over Correlated 
Not-Necessarily Identically Distributed Fading Channels 

We utilize the MacRobert's E function representation of the UP expression (i.e., see Corollary [3]) 
instead of the other representations. More specifically, using the well-known the integral formula of the 
MacRobert's E function d Eq. (6)], §\ Eq. (7.814/2)], given by 



exp (-n) E {a,}U {p^Y,^,; - du = E [{a,}th {pAU ^] 



(C.l) 



where 3ft {ttp+i} > and z E M^, and performing some algebraic manipulations, one can readily obtain 
an alternative representation of (IC.ll) as 



(«7, 



end) 



E 



du. 



(C.l) 



Accordingly, substituting (|C.2I) and the total instantaneous SNR 7g„^ = ^^^^ 7^ into it is straightfor- 
ward to show that the AUP expression Paup = '^[Pup (lend)] L-branch diversity combiner can be 
explicitly given by 



n 



Paup = 1 " 3 ^ ^ 



m 



exp 



u 



aEf=i7£, 



E 



A 



6+1; 



u 



du 



(C.3) 



Note that, E 



exp(— s/^^^]^ 7^) with s = u/a can be considered as the MGF of the reciprocal distribution 
lend = J2e=i 7^- order to proceed further, we present the following lemma. 

Lemma 1 (MGF of Reciprocal Distribution [|20l Theorem 1]). Given any nonnegative continuous RV TZ 
distributed over (0, 00) with the PDF pnij") <^f^d the MGF A^7^(s) for 3ft{s} > 0. Then, the MGF of its 
reciprocal distribution TZ (i.e., TZ = 1/71) is given by 

d 



du 



Mniu) 



du 



(C.4) 



with the convergence region 3ft{s} > 0, where Jq (•) is the zeroth-order Bessel function of the first kind 
defined in [[6l Eq. (6.19.7)]. 
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Accordingly, applying Lemma[I]on the expectation part E exp(— (n/a)/^^^^^ 7^) of (IC.3I ) results as 



n poo foo 




u^-\j^i2j^ 1 X 




E 



u 



du 



|exp^-pX:7. 



dp). (C.5) 



Note that using the Hankel transform of the MacRobert's E function [F, Eq. (7.823/1)], and then applying 
ETl Eq. (1)], the inner integral of (IC.5I) can be easily obtained in closed-form as 



du = G^'"i 2 



6,0 



(C.6) 



Finally, replacing the inner integral in (IC.5I) with the Hankel transform of MacRobert's E function given 
by (IC.6I) and using the definition of joint MGF (i.e., M.j^^^^^,„^j^ (p) = E[exp(— p 7^)]), (IC.5I) simplifies 
to (|301 ), which proves Theorem [3l 
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TABLE I 

Parameters a, b, n and d for the BEP and Capacity Performance Measures 



Performance Measure 


a 


b 


n 


d 


BEP of orthogonal coherent BFSK (T] Eq. (8.43)], 

Pbep (lend) ~ Q iV^Und) ' where Q (•) is the Gaussian Q-function having a one-to-one mapping 
with the complementary error function erfc (■), i.e., Q (2:) = ^erfc(.2:/\/^). 


1 
2 


1 

2 


1 


00 


BEP of orthogonal noncoherent BFSK Q] Eq. (8.69)], 

PBEp{l^„d) = iexpf-^y 


1 
2 


1 


1 


00 


BEP of antipodal coherent BPSK [JJ Eq. (8.19)], 

i^BEP (7end) " "-i \V /end) ■ 


1 


1 

9 


1 


00 


BEP of antipodal differentially coherent BPSK (DPSK) (T| Eq. (8.85)], 

Pbep {lend) = ^ cxp (-7e„d)- 


1 


1 


1 


00 


BEP of correlated coherent binary signaling fTl Chapter 8, footnote 6], 

Pbep (lend) = Q{V^<^%nd)- 


a e [0,1] 


1 

2 


1 


CO 


Shannon Capacity [T] Eq. (15.22)], 

Pc (lend) ~ '°§ (1 + '^lend)' where a is the transmitted power and log (■) is the natural logarithm 
(i.e., the logarithm to the base e) (6| Eq. (4.1.1)]. 




1 


2 


1 



TABLE II 

MGFs OF SOME Well-known Fading Channel Models 



Instantaneous SNR Distribution, i.e., (7) 



MGF M^^ (s) and its derivative —M^^ (s 



One-Sided Gaussian fH Sec. 2.2.1.4] 



1 



2V^7, 



■H, 



0,1 



" 7 









where 7^ is the average power (i.e., 7^ > 0). 



M,, is) = ^H};} 



1 


(1,1)" 


_27^s 


(ii). 



1 

v/TT 



1 


1 ■ 


1 


_27^s 


1 

2 - 


^1 + 27,s 



ds 



-M,^ is) = ^H^;^ 



27^5 



(2,1), (1,1) 
(1,1), (2,1) 



2,1 



1 


2,1" 


_27^s 


2' ^ J 



7£ 



(1 + 27fS) 



3/2 ' 



where G™q"[-] and H™g"[-] represent the Meijer's G function Q EQ- (8.2.1/1)] and Fox's H function Eq. (8.3.1/1)], respectively. 
Note that one-sided Gaussian fading coincides with the worst-case fading or equivalently, the largest amount of fading (AoF) for all 
Gaussian-based fading distributions. 



Exponential H] Eq. (2.7)] 



Vi, 7 = — cxp ' 

In 



-H 



It 



0,1 



where 7^ is the average power (i.e., 7^ > 0). 



(0,1) 



-^7, is) - H 



1,1 



(1,1) 
(1,1) 



'-'1,1 



'7£ 



^-^7.(-^)=7£H^;2 



'7^ 



(2,2), (1,1) 
(2,1), (2,1) 



2,2 



1 



2,1 
2,2 



It 



(1 + ^7,) 



Gamma HI Eq. (2.21)] 

Pi, (7) 



r(mj V 7£ 



exp 



TO^7 



r(m^)7^ 



ttI.O 



L7f 



^7 



(m^ - 1,1) 



where 7^ is the average power, and where (0.5 < m^) 
denotes the fading figure. Moreover, r(-) is the Gamma 
function (8) Sec. 8.31]. 



-^7, (■') 



r(m^) 



-H 





(1,1) ■ 


^1,1 


rrtj! 


1 


7^s 


(m^, 1) 


r(m^ 


) 



7f 



H 



2,2 



(2,1), (1,1) 
K + 1,1), (2,1) 



G 



2,1 





2,1 




m^-l-1,2 



r(m^ + 1) 



-7.(1 + -^ 



Note that the Nakagami-m distribution spans via the m parameter the widest range of amount of fading (AoF) among all the multipath 
distributions [T]. As such, Nakagami-(j' (Hoyt) and Nakagami-n (Rice) can also be closely approximated by Nakagami-Hi distribution [T] 
Eq. (2.25)], (U Eq. (2.26)]. 



WeibuU d Eq. (2.27)] 



cxp 



It 



(i-i/e„i/e,) 



where = r(l + l/5£) and where (0 < 5^) denotes the 
fading shaping factor. Moreover, 7^ is the average power. 



-^7, (5) = H};} 



(1,1) 
(l,f 



2;,2fe 

2fc+2/-2 ^2fc,2/ 



— TW^, s) = H2 2 




2k 



^(-2fc) 
'"(2fc) 

"(1) 
"(20 



VieFI 



(l + i,i),(2,l) 



n2/£+2;-2 
(iTT) S 



-G: 



2;-i-i,2fe 

2fe-|-l,2J-|-l 



c-jf (2fc) 



2fc 



s^f'Vil'i{2l)pl 



'"(2fc) 
"(1) 1 



where the Meijer's G representations are given for the rational value of the fading shaping factor (that is, we let = k/l, where fc, 
and I are arbitrary positive integers.) through the medium of algebraic manipulations utilizing Q Eq. (8.3.2.22)]. In addition, the coefficient 

i±l, . . . , 5^±2i^ with a; e C and n e N. 



(x) (x) 

H) . of the Meijer's G function is a set of coefficients such that it is defined as H; . 

in) J in) 



TABLE III 

MGFs OF SOME Well-known Fading Channel Models 



Instantaneous SNR Distribution, i.e., (7) 



MGF (s) and its derivative —M-y, (s 



ds 



Hyper-Gamma 

K 

Pit (7) 



Eq. (3)] 



^1 r(TO^^,) V7£fc 



i-ek 



e >fc , 



K 

E 

k=l 



0,1 



"•-Ik 
Ilk 



7 



(m 



1,1) 



where m^^, (0.5 < rn^^) is the fading figure, 7^^ (0 < 7^^) 
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the zeroth order modified Bessel function of the first kind (6] 
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where ni (0 < n^) and 7^ (0 < 7^) are the LOS figure and 
average power, respectively. 
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fc + 1 and f2j. = (fc + 1) -^^^^a , respectively. In addition, the weighting coefficients are given 

by = '^i' Gxp (— ra^) /r(A; + 1). It may be useful to notice that the line-of-sight (LOS) figure i.e. rig is related to the Rician Ki 
factor by Kg = n| which con'esponds to the ratio of the power of the LOS (specular) component to the average power of the scattered 
component. 
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where /^^(dB) and a'^(dB) are the mean and the standard 
deviation of 7„. 
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(abscissas) of the X-order Hermite polynomial (6] Table 25.10]. 
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where m^^ < nJ^^) denotes the shadowing severity, and 7^^ (0 < 7^^) represents the average power. In addition, (■) is the nth 
order modified Bessel function of the second kind (6) Eq. (9.6.24)]. 
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where <j)t = f^^^ '^''^e V'^ = '^sf ~ ^nd where the parameters (0.5 < m^) and m^^ (0.5 < ti^^) represent the fading figure 
(diversity severity / order) and the shadowing severity, respectively. 7^^ (0 < 7^^) represents the average power. 
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where the parameters (0.5 < < oo) and (0 < < oo) represent the fading figure (diversity severity / order) and the 
fading shaping factor, respectively, while 7^ (0 < 7^ < oo) is the average power In addition, the parameter is defined as = 
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where <S>i = 2-K{lks)™-t (Zsfc)'"sf/V(27r)'*=+'s'=+*^'==-iZ«s and 'I'f = (ifcs)'^* {lsk)'-st /{kks). and where the parameters (0.5 < 
< 00) and (0 < < 00) represent the fading figure (diversity severity / order) and the fading shaping factor, respectively, while 
(0.5 < m^^ < 00) and (0 < < 00) represent the shadowing severity and the shadowing shaping factor (inhomogeneity), 
respectively. In addition, the parameters and f}. are defined as = F (m^ + 1/5^) /F (m^) and = F {m^p + l/^sf) ('"sf)' 
respectively, where F (■) is the Gamma function (6] Eq. (6.5.3)]. In addition, F (■, ■, ■, ■) is the extended incomplete Gamma function defined 
as F (a, X, b, 13) = exp (— r — br~/^^ dr, where a, /3, b 6 C and x S M"*" M7\ Eq. (6.2)]. Referring the coefficients of the 

Meijer's G function, sj^? is a set of coefficients such that it is defined as s!^? = — , 2LLi , , x+n—i ^j^j^ x £ C and n S N. 
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Fig. 1. ABEP performance comparison of BPSK and BDPSK binary modulations over generalized Gamma fading channel with the 
parameters: (a) = 2, = 0.25, (b) = 2, = 0.5, (c) = 2, ^£ = 0.75, (d) rrig = 2, = 1.0 and (e) = 2, = 1.5. The 
number of samples is chosen as 10^ in the computer-based simulations. 
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Fig. 2. ABEP performance comparison of BFSK and NC-FSK binary modulations over generalized Gamma fading channels with the 
parameters: (a) = 2, = 0.25, (b) = 2, = 0.5, (c) = 2, = 0.75, (d) = 2, = 1.0 and (e) = 2, = 1.5. The 
number of samples is chosen as 10^ in the computer-based simulations. 



IEEE TRANSACTIONS ON COMMUNICATIONS, VOL. X, NO. XX, JAN. 2011 



28 




Fig. 3. AC of wireless communication system over generalized Gamma fading channel with the parameters: (a) — 2, — 0.25, (b) 
rrii = 2, ^£ = 0.5, (c) = 2, = 0.75, (d) = 2, = 1.0 and (e) = 2, = 1.5. The number of samples is chosen as 10*' in 
the computer-based simulations. 
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Fig. 4. ABEP performance comparison of BPSK and BDPSK binary modulations for L-branch MRC receiver over generalized Gamma 
fading channels with the parameters — 2 and = 0.25. The number of samples is chosen as 10^ in the computer-based simulations. 
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Fig. 5. ABEP performance comparison of BFSK and NC-FSK binary modulations for L-branch MRC receiver over generalized Gamma 
fading channels with the parameters = 2 and — 0.25. The number of samples is chosen as 10' in the computer-based simulations. 
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Fig. 6. AC of wireless communication systems witli an L-brancli MRC receiver over generalized Gamma fading channel with the parameters 
TTif = 2, = 0.25. The number of samples is chosen as 10*^ in the computer-based simulations. 



